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' >' Abstract 

Ph Recent developments on financial markets have revealed the limits 

^S^ of Brownian motion pricing models when they are applied to actual 

d markets. Levy processes, that admit jumps over time, have been found 

H more useful for applications. Thus, we suggest a Levy model based 

'— ' on Forward-Backward Stochastic Differential Equations (FBSDEs) for 

p^ option pricing in a Levy- type market. We show the existence and 

J> uniqueness of a solution to FBSDEs driven by a Levy process. This 

^V result is important from the mathematical point of view, and also, 

1/-^ provides a much more realistic approach to option pricing. 
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1 Introduction 

Since the seminal contribution made by F. Black and M. Scholes [2], several 
methodologies to value contingent assets have been developed. From Plain 
Vanilla options to complex instruments such as Collateral Debt Obligations 
or Baskets of Credit Default Swaps, there are models to price virtually any 
type of contingent asset. There are, indeed, successful attempts on providing 



general models which in theory could price any kind of contingent claim 
(see, for example [3]), given a payoff function. The idea behind these models 
is quite standard: A portfolio replicating the payoff function of the asset 
is devised and, under non-arbitrage conditions, the price of the asset at a 
certain instant of time is the price of the replicating portfolio at this time. 

However, in spite of all this diversity and sophistication, there is an as- 
sumption that is, up to recent times, rarely questioned. Specifically, we refer 
to the assumption that stock prices are continuous diffusion processes, pre- 
supposing thereby that the returns have normal distributions at any time. 
However, it is well known today that empirical distributions of stock prices 
returns tend to deviate from normal distributions, either due to skewness, 
kurtosis or even the existence of discontinuities (Eberlein et al. give evidence 
of this phenomenae in [4]). The recent developments have shown that the 
reliance on normal distribution can bring costly surprises, especially when 
extreme and disruptive events occur with a much higher frequency than that 
estimated by models. 

As such, we believe that no matter which historical status the normal 
distribution has acquired throughout the years, strong efforts should be un- 
dertaken in order to develop alternative models that incorporate assumptions 
adequated to the observed evidence on financial markets, such as asymme- 
try or skewness. We do not pretend that some definitive model can actually 
be developed, especially when market participants' main activities are cur- 
rently shifting due to the conditions imposed on financial markets. Indeed, 
the recently introduced regulations on financial markets severely restraining 
the use of own's capital for trading purposes will force the financial players 
to find new ways of driving a profit. This adds another layer of uncertainty 
about the assumptions imposed on a model. We believe, however, that in 
spite of the inherent inability to prove that any present model can account 
for future market conditions, it is worth to attempt to correctly price finan- 
cial claims in the present and near future market conditions, which, as it is 
clear now, is fundamental to the stability of markets. 

Taking the above considerations into account, we propose to replicate 
contingent claims in Levy-type markets, i.e. in markets with the stock-price 
dynamics described as St = Soe^\ where Xt is a Levy- type stochastic in- 
tegral defined in [1]. This allows the relaxation of conditions posed on the 
pricing process such as symmetry, non-skewness or continuity, imposed by 
the Brownian framework. The self-similiarity of the pricing process, appear- 
ing due to a Brownian motion, is also ruled out from the assumptions. We 



base our model on the study of Forward-Backward Stochastic Differential 
Equations (FBSDEs) driven by a Levy process. FBSDEs combine equations 
with the initial and final conditions which allows one to search for a replicat- 
ing portfolio. Specifically, we are concerned with the following fully coupled 
FBSDEs: 
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* ^ Yo V V \ ^u(^) 



+ /; fis, X„ Ys, Zs) ds + YZi /o ^^(^' ^-' Ys-) dm 



h{XT) + j^ g{s, X„ Y„ Zs) ds - YZi !l Zl dHf\ ^^'^^ 

where the stochastic integrals are written with respect to the orthogonalized 
Teugels martingales {Hi }^^ associated with a Levy process Lt [11]. We are 
searching for an M^ x M9 x {MP x £2)-valued solution (Xj, Yt^ Zt) on an arbi- 
trary time interval [0, T], which is square-integrable and adapted with respect 
to the filtration J^t generated by Lt- To the authors' knowledge, fully coupled 
FBSDEs of this type have not been studied before. Fully decoupled FBSDEs 
involving Levy processes as drivers were studied by Otmani [6]. Backward 
SDEs driven by Teugels martingales were studied by Nualart and Schoutens 
[10] . Our method of solution to the FBSDEs could be compared to the Four 
Step Scheme [7]. The original four step scheme deals with FBSDEs driven 
by a Brownian motion, and the solution is obtained via the solution to a 
quasilinear PDE. Replacing the stochastic integral with respect to a Brow- 
nian motion with a stochastic integral with respect to the orthogonalized 
Teugels martingales leads to a partial integro-differential equation (PIDE). 
The solution to the PIDE is then used to obtain the solution to the FBSDEs. 
The organization of the paper is as follows. In Section 2, we give some 
preliminaires on the martingales {Hi }. In Section 3, under certain assump- 
tions, we prove the existence and uniqueness result for the associated PIDE. 
Our main result is Theorem 3.7, where we obtain a solution to FBSDEs (1.1) 
via the solution to the PIDE and prove its uniqueness. In section 4, we ap- 
ply the results of Section 3 to model hedging options for a large investor in 
a Levy-type market. Previously, this problem was studied by Cvitanic and 
Ma [3] for a Brownian market model. Finally, we study conditions for the 
existence of replicating portfolios. 

2 Preliminaires 

Let {Q,J-',J-'t,P) be a filtered probability space, where {J^t}, t G [0,T], is 
the filtration generated by a real- valued Levy process Lt. Note that the Levy 



measure v of Lt always satisfies the condition 

(1 Nx^')v{dx) < oo. 



We make the filtration J-'t P-augmented, i.e. we add all P-nuU sets of J-" to 
each J^f. Following Nualart and Schoutens [9] we introduce the orthogonalized 
Teugels martingales {H^}°Zi associated with Lt. For this we assume that 
for every e > 0, there exists a A > so that 

exp(A|x|) z/((ix) < oo. 

The latter assumption guaranties that 

|a;|* ^{dx) < oo for z = 2, 3, . . . . 



It was shown in [11] that under the above assumptions one can introduce 
the power jump processes and the related Teugels maringales. Futhermore, 
it was shown that the strong orthogonalization procedure can be applied to 
the Teugels martinagles and that the orthonormalization of the Teugels mar- 
tingales corresponds to the orthonormalization of the polynomials 1, x, x^, . . . 
with respect to the measure x'^v{dx) + a^(5o((ix), where the parameter a G M 
is defined in Lemma 2.1 below. As in [10], by {qi{x)} we denote the system 
of orthonormalized polynomials such that gj_i(x) corresponds to H). . Also, 
we define the polynomial Pi{x) = xgj_i(x). We refer to [11] for details on the 
Teugels martingales and their orthogonalization procedure. In the following. 
Lemma 2.1 below will be usefuU. 

Lemma 2.1. The process H^ can be represented as follows: 



i/f = g,_i(0)P\t) + [ p,{x)N{t,dx), 

Jr 



where B^(t) = ^j=iAj-Bj(t) with X^\ = a, \i E M, {Bi{t)}f^^ are indepen- 
dent real-valued Brownian motions, and N{t, A) is the compensated Poisson 
random measure that corresponds to the Poisson point process ALt. 

Proof. We will use the representation below for H^ obtained in [10]: 



iff = g._i(0)L, + J2 M^^s) -tE Y, p{AL 



0<ss£t 



0<ss£l 



-tg,_i(0)E[L 



ij' 



where Pi{x) = Pi{x) — xgj_i(0), and E is the expectation with respect to P. 
Taking into account that Lt = L^ + J2o<s<t'^^s, where L^ is the continuous 
part of Lt, we obtain: 



^f = q^^l{o)^ + Yl p*(^^^) -^^ Yl p'(^^ 



0<s^t 



0<s^l 



%_i(0)S[Li 



g,_i(o) \li - Em]] + J2 P^i^^s) -E Yl Pi(^L^) 



0<ss£t 



0<ss£i 



= qi_,(0)B\t) + / pi{x)N{t,dx). 

Jr 

In the sequence, the following lemma will be frequently applied: 
Lemma 2.2. It holds that 

Pi{x)pj{x)u{dx) = Sij - a^gi_i(0)gj_i(0). 
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Proof. The proof is a straightforward corollary of the orthonormality of 
qi-i{x) with respect to the measure x'^u^dx) + a^6o{dx). D 

We will need an analog of Lemma 5 from [10] which was proved in the 
latter article for a pure-jump Lj. We obtain this result for the case when Lt 
has both the continuous and the pure-jump parts. 

Lemma 2.3. Let /i : fi x [0, T] x M — > M" be a random function satisfying 



E [ \h{s,y)\'iy{dy) <oo. 
Jo 



(2.1) 



Then, for each t G [0, T], 



Y HsALs) = Y f I u{dy)h{s,y)p,{x)dHf+ I j h{s,yMdy)ds. 



Proof. Note that 

Mt= y^ K^^^Ls)- / / h{s,y)u{dy)ds = / / h{s,x) N{ds,dx) 
Q<^<j Jo Jr Jo Jr 

(2.2) 

is a square integrable martingal, i.e. sup^g^.T] -E'|Mip < oo, by (2.1). By the 
predictable representation theorem [9], there exist predictable processes ipi 

/o Zli^i I'^i? < ^ aiid such that Mt = J2'iLi /o fi{s)dHs . Since 



with E 



/o 



{H^'\H^^^)^ = t6ij [11], then 

(M,/7«)^= / ip,{s)ds. 
Jo 

On the other hand, by (2.2) and Lemma 2.1, 



{M,H^'^) 



h{s,x)N{ds,dx),qi-i{0)Bt+ / pi{x)N{dt,dx) 
Jr Jo Jr it 

/ h{s,x)pi{x)iy{dx)ds. 

Jr 

This implies that (pi{s) = J^h{s, y)pi{y)h'{dy) , and therefore, 

V/i(s,AL,)-/ [ h{s,yMdy)ds = f2 [ [ h{s,y)p,{yMdy)dHi^^ 
Ks^t -^0 Jr ^^1 Jo Jr 



0^ 
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3 FBSDEs and the associated PIDE 

3.1 Problem Formulation and Assumptions 

Consider the FBSDEs: 



Xt = x + /J f{s, X„ n, Zs) ds + j; ais, Xs-, n_) dH,, 

Yt = HXt) + J^gis,Xs,Ys,Z,) ds - J^ Z, dH,, (3.1) 

tG[0,T], 



where 



[0, T]xR^ xR^ X (M^ X £2) -^ M^, 
[0, T]xR^ xR^ ^ 



)P 



[o,r] X M' X M-* X 



)Q 



X£2)^ 



— )■ 



)Q 



are Borel-measurable functions. Here, the stochastic integrals 

/ a{s,Xs-,Ys-)dHs and / ZgdHs 
Jo Jt 



are shorthand notation for 

00 ^t 



Y, a,{s,X,_,Y,_)dHf^ and ^/ Z^ rfi7» 
i=i -^0 i=i Jt 

respectively, where Z, = {ZI}°1^, a = MZi, (^i ■ [0,T] x R-^ x M^ ^ R^. 
The solution to FBSDEs (3.1), when exists, will be an M^ x R^ x (M^ x £2)- 
valued J-j-adapted triple {Xt, Yt, Zt) satisfying 



•^0 ^ i=i 



(it < 00, 



and verifying (3.1) P-a.s.. The latter includes the existence of the stochastic 
integrals in (3.1). Implicitly, we are assuming that Xt and Yt have left limits, 
and that Zt is J-t-predictable. So in fact, we are searching for cadlag (X^, Yt), 
which will guarantee the existence of Xt- and Yt-, and predictable Zt. 
We associate to (3.1) the following final value problem for a PIDE: 

{dtd{t, x) + f\t, X, 9{t, x), d^^\t, x)) dkOit, x) + I3^\t, X, 9{t, x)) dl 9{t, x) 
- 4 [e{t, X + 6{t, X, e{t, x),y)) - e{t, x) - duOit, x)5\t, X, e{t, x), y)] v{dy) 

+g{t,x,e{t,x),e^^\t,x)) = 0, 

{e{T,x) = h{x) 



(3.2) 



with ^(1) : [0, T] 



PQ 



xh, 



ef\t,x)= / [e{t,x + 5{t,x,9{t,x),y)) - e{t,x)]pi{y)v{dy) 



+ c'^{t,x,e{t,x))dke{t,x). (3.3) 



The connection between /3^', 6, cf and the coefficents of FBSDEs (3.1) is the 
following: 



S{t,x,y,y') = ^ai{t,x,y)pi{y'), 



(3.4) 



a 



2 °° 



/3'='(t,x,i/) = y(^^af(t,x,y)g,_i(0))(5^a5.(t,x,i/)g,_i(0)), (3.5) 
cf(t,a;,y) = (Tf(t,x,?/) - / 5''{t,x,y,y')pi{y')u{dy'). 



(3.6) 



To guarantee the existence of the above functions we will make the assump- 
tion: 

Ao J:T=l<l^-M'<^■ 

Since a^ = {af}°Zi takes values in £2, AO immediately guarantees the conver- 
gence of the both multipliers in (3.5). The convergence of the series in (3.4) 
is understood in L2{i'{dy')) for each fixed (t,x,y). Moreover, it holds that 



2 

^ai{t,x,y)pi{y') v{dy' 



1=1 



|cr(t, X "'^"^ 



'W/'llR^x^a 



-a^^ai(t,x,|/)gj_i(0) . (3.7) 



j=i 



Indeed, applyng Lemma 2.2, for each fixed A^ we obtain: 

I \^(^i{t^^^y)Pi{y') ^{dy') = X^(c^i,^i) / Pi{y')Pj{y')^{dy') 

N N 



i=l 



Now letting A^ tend to infinity, we obtain (3.7). 
Lemma 3.1. The following assertions hold: 

1. c'l{t,x,y) =a'^qi-i{0)J2T=iCrj{s,x,y)qj^i{0). 



2. For each k, c'' = {c^}°Zi takes values in £2- 

3. For each {s,x,y), \ f^6{s,x,y,y')pi{y')u{dy')\ e £2- 



00 



Proof. Define 6]\f{s^x,y,y') = ^j=iCrjis,x^y)pj{y'). By what was proved, 
for each {s,x,y), 6iy{s,x,y, ■) — )■ S{s,x,y, ■) in L2{i'{dy')), and therefore, 
for each i, and for each (s, x, y), 



6N{s,x,y,y')pi{y')u{dy') ^ / 6{s,x,y,y')pi{y')u{dy') 

Jr 

as A^ — > 00. On the other hand, by Lemma 2.2, 

/ SN{s,x,y,y')pi{y')u{dy') = ai{s,x,y) - a'^ qi_i{0)^aj{s,x,y)qj^i{0). 

J* ,=1 

Comparing the last two relations, we obtain that 

« 00 

/ S{s,x,y,y')pi{y')u{dy') = ai{s,x,y) - a^ qi_i{0)^aj{s,x,y)qj_i{0) 

Jr ^.=1 

(3.8) 

which proves Assertion 1. Assertion 2 is implied by Assumption AO and 
Assertion 1. Finally, (3.6) implies Assertion 3. D 

The heuristic argument behind PIDE (3.2) assumes the connection Yf = 
6{t,Xt) between the solution processes Xt and Yj to (3.1) via a C^'^-function 
6. Ito's formula applied to 6{t,Xt) at points t and T leads to another BSDE 
which has to be the same as the given BSDE in (3.1). Thus we "guess" PIDE 
(3.2) by equating the drift and stochastic terms of these two BSDEs. 

3.2 Solvability of the PIDE 

We solve Problem (3.2) for a particular case when function f{t,x,y,z) 
does not depend on z, and g(t,x,y,z) = gi(t,x,y) + {g2it,x,y),z), where 
92{t,x,y) = {g2\t,x,y)}^^ G £2 for each {t,x,y), and {g2{t,x,y),z) = 

X^i^i 92 v-1 ^-1 y) ^i- Thus, we are dealing with the following final value prob- 
lem for a PIDE: 

dtOit, x) = -[A{t, e{t, ■))e]{x) + gi{t, X, e{t, x)), 

6{T,x) = h{x), 



9 



where A{t,p{t, ■)) is a partial integro-differential operator given by 

[A{t,p{t,.))e]{x) 

= [f + {92, c"),,] (t, X, p{t, x)) dke{t, x) + I3^\t, X, p(t, x)) dl e{t, x) 

(3.10a) 

+ / [9(t,x + 5{t,x,p{t,x),y)) - 9{s,x) - dk9{t, x)5^{t, x, p{t, x),y)] u{dy) 
Jr 

(3.10b) 

CO „ 

-X^5'?^(^>a;,p(t,a;)) / [9{t,x + S{t,x, p{t,x),y)) - 9{t,x)]pi{y)u{dy) 

(3.10c) 

with the domain D{A{t, p{t,-)) = C^(]R^ -^ R^), the space of bounded 
continuous functions M.^ — )■ M*^ whose first and second order derivatives are 
also bounded. We assume the following: 

Al Functions /, gi, g2, and a are bounded and have bounded spatial deriva- 
tives. 

A2 One of the following is fulfilled: 

• ^2 = 0, 

or 

• lR\y\^{dy) < 00, T.Zi\92\t,x,y)\J^\pi{y')\iy{dy') < Kt, and 

1 

j:T=i\9^\t,^,y)\{L\My')\Mdy')y < K^, where K,<oo de- 
pends only on t. 

Lemma 3.2. Let Assumptions AO, Al, and A2 be fulfilled. Then A{t, p{t, ■)), 
defined by (3.10), is a generator of a strongly continuous semigroup on 
Cb{R^ ^R^). 

Proof. Note that by Assertion 1 of Lemma 3.1 and by (3.5), functions c'^ and 
/3'^' are bounded and Lipschitz in the spatial variables. This implies that the 
SDE 

00 
dX^ = [f + (^2, c%,) {t, X„ p(t, Xt)) + Y, ^Ht, Xt., p{t, Xt.))dHf\ 

i=l 

(3.11) 



10 



with Xq = X, has a pathwise unique cadlag adapted solution. The existence 
and uniqueness of a solution to an SDE of type (3.11) will be proved in 
Paragraph 3.3. Now application of Ito's formula to ip{Xt), where ip is twice 
continuously differetiable, shows that the operator given as the sum (3.10a) 
+ (3.10b) is the generator of the solution to SDE (3.11), and therefore, it 
generates a strongly continuous semigroup on Cb(]R^ — )■ MP). Now ii g2 = 0, 
then (3.10c) is zero. Otherwise, as a function of x, it belongs to the class 
CbiR^ -^ R^). Indeed, assume first that 9{t, ■ ) e C^(M^ ^ M«), and prove 
that the map 



dP 



^R'^, x^ e{t,x + 6{t,x,p{t,x),y))pi{yMdy) (3.12) 



is Lipschitz, where x and t are fixed. Fix xi and X2 from M^, and let pi 
p(t,xi) and p2 = p(t,X2), where t is fixed. We have: 



[9{t, X + 6{t, xi, pi, y) - 9{t, X + 6{t, X2, P2, y)]pi{y)v{dy) 

^ max|Ve(t,x)| / \5{t,xi,pi,y)-5{t,X2,p2,y)\\pi{y)\y{dy) 
^ max|V6'(t,a:)|( / \5{t,xi, pi,y) - 5{t,X2, p2,y)?T^{.dy) / \p,{y)\^v{dy) 



^ Kmax\Ve{t,x)\[ / \pi{y)fjy{dy)) \\a{t,xi,pi) - a{t,X2,p2)\\RPxe2^ 

(3.13) 

where i^ is a constant determined by (3.7), and is the same for all i. The 
Lipschitzness of map (3.12) for 6{t, ■ ) G C^(]R''^ — t- M*^) follows from the Lips- 
chitzness of the map a(t, ■ , ■ ). The boundedness and the continuity of (3.10c) 
are implied by Al and A2 and by the continuity of 9(t, x) and g2(t, x, p{t, x)) 
w.r.t. X. If 9(t, ■ ) does not possess a continuous bounded derivative, then we 
can find a.9 E C^(M^ — >■ M*^) with sup^jg^p \9(t,x)—9{x)\ < e, for an arbitrary 
small e > 0. By Assumption A2, (3.10c) will be of the class Ch{R^ -^ M^) 
for each fixed t. Now A(t,p(t, ■ )), defined by (3.10), is a generator as a sum 
of a generator and of a bounded operator on Cb(]R^ — > MP). D 

The common method to deal with problems of type (3.9) is to fix a C^' - 
function p{t, x), and consider the following non-autonomous evolution equa- 
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tion: 

dtd{t,x) = -[A(t,p(t,-))^](x) - (7i(t,x,p(t,x)), 

9{T,x) = h{x). ^ ■ ' 

Suppose h G C^(]R'^ — t- IR*^), i.e. it belongs to the common domain of the 
generators A{t, pit,-)). Then, there exists a backward propagator U{s,t,p) 
[5], ^ s ^ t ^ T, so that 

e{t,x) = [U{t,T,p)h]{x) + J [Uit,s,p)g,is, ■,p{s,-))]ix)ds. 

We organize the map 

$ : Cfe([0, T] X M^ ^ R^) -^ Cfe([0, T] x M'f' ^ M^^)^ ^ ^ q^ (3_X5) 

and prove the existence of a fixed point. We would like to use the following 
result from [8]. 

Lemma 3.3 (Convergence of backward propagators, Kolokoltsov). Let E he 
a Banach space. Suppose we are given a sequence of backward propagators 
U"'{s,t) generated by the families A^ and a backward propagator U{s,t) gen- 
erated by the family At, and that the generators A^ and At have the same 
common domain D. Further we suppose that all these propagators satisfy the 
conditions formulated in Theorem 1.9.4 pO'g^ 55 of [8] and that all U{s, t) are 
uniformly bounded as operators in D. Then for any (f & D, 

sup ||f/"(t,s)v9-t/(t,s)v5|U ^ K ess supt^s^rW^"^ - ^s\\c{d,e)\\^\\d- 

(3.16) 

Remark. In [8] inequality (3.16) is stated in the form 

sup W^it, s)ip - U{t, s)ip\\E ^ K ess supo^j.^^ \\At - At\\c(D,E)\\v\\D- 

(3.17) 

However, analysis of the proof of the above statement in [8] shows that (3.16) 
can be obtain in exactly the same way as (3.17). 

Define E = Cb(R^ -^ R^) and D = ClCR^ -^ R^). 
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Lemma 3.4. Let Assumptions AO, Al, and A2 hold. Then, for any function 
if E E, there exists a constant K > so that 

sup \\U{t,s,p)ip-U(t,s,p)ip\\E ^ K sup \\p{s,x) - p{s,x)\\e\\(p\\e- 

s£[t,T] s£[t,T] 

Proof. It is known that ess sup^ ||/(a^)|| ^ sup^ 11/(2^)11 ^^i any Banach space- 
valued function /. Therefore, 

ess supt^,^-rp(s,p(s, ■)) - A{s, p'{s,-))\\c{d,e) 

^ sup \\A{s,p{s,-))-A{s,p'{s,-))\\ciD,E)- (3.18) 

We would like to prove that the norm on the right-hand side is smaller than 
\\p{s, ■ ) — p'{s, ■ )\\e- Taking into account that 

\\e{s, ■)\\d = sup |^(s,x)| + sup |V^(s,x)| + sup |VV^(s,x)|, 



and applying (3.10), (3.5), and Lemma 3.1, we obtain that there exists a 
constant K > which does not depend on s, so that 

sup sup \A{s, p{s,x))9 — A{s, p'{s,x))6\\ 



Id^Ix 



^ K sup 



\f{s,x,p{s,x)) - f{s,x,p'{s,x))\ 

+ \\92is,x,p{s,x)) - 5-2(5, x,p'(s,x))||£2 

+ \Wis,x,p{s,x)) - o-(s,a;,p'(s,a;))||]RPx£2 

+ sup \e{t,x')\\\g2is,x,p{s,x)) - g2is,x,p{s,x))\\i^ 



+ sup |V^(t,x')| y^\gt\s,x,p{s,x))\{ \ \v^{y)fy{dy) 

r. 1 

X ( / |5(t,a;,p(t,x),y)-5(t,x,p'(t,a;),y)|V(dy))' (3.19a) 



x<& 



+ sup |VV^(t,x')| / \b{t,x,p{t,x),y) - b{t,x,(I{t,x),y)fv{dy) 



(3.19b) 



To obtain (3.19a), we applied the same argument as in (3.13) to the differ- 
ence of the terms of type (3.10c), and to obtain (3.19b), we applied Taylor's 
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formula to the difference of the terms of type (3.10b). (3.13) also shows that 
(3.19a) is smaller than \\a{s,x, p{s,x)) — cr(s,a;,p'(s,x))||iRPx^2 ^P ^^ ^ mul- 
tiplicative constant. The same applies to (3.19b) due to (3.7). Therefore, 
modifying the constant K, if necessary, by Assumption Al, we obtain that 

sup \\A{s,p{s, ■ ))9 - A{s,p'{s, ■ ))9\\e ^ K\\p{s, ■ ) - p\s, ■ )IU 

where the constant K does not depend on s. Now by (3.18) and Lemma 3.3, 

sup \\U{t, s, p)ip - U{t, s, p')ip\\E ^ K snp ||p(s, ■ ) - p'(-5, ■ )IU llv^lU 

s&[t,T] s&[t,T] 

which completes the proof. D 

Theorem 1. Let Assumptions AO, Al, and A2 hold. Then, the problem (3.9) 
has a unique mild solution, i.e. there exists a unique function 9 G Cb([0, T] x 
M^ ^ M^), so that 

9{t,x) = [U{t,T,9)h]{x)+ f [U{t,s,9)gi{s,-,9{s,-))]{x)ds. (3.20) 

Proof. The proof of the existence and uniqueness of a solution to (3.20) is 
equivalent to the existence of a unique fixed point of map (3.15). The latter 
follows from Lemma 3.4 and from Assumption Al relative gi. D 

Corollary 1. Suppose AO, Al, and A2 hold, h e C^iR^), and the solution 
9{t, ■ ) to (3.20) belongs to Cg(M^ -^ R^) for each t G [0,T]. Then, 9 is also 
a solution to (3.9). 

3.3 Existence and Uniqueness Theorem for the 
FBSDEs 

In Paragraph 3.2 we found some conditions under which there exists a unique 
solution to PIDE (3.2). However, this solution may exist under more general 
assumptions. Thus, we prove the existence and uniqueness of a solution to 
FBSDEs (3.1) assuming the existence and uniqueness of a solution to PIDE 
(3.2). We will assume the following: 

A3 Functions /, g, and a possess bounded first order derivatives in all 
spatial variables. 
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A4 Assumption AO is fulfilled and Final value problem (3.2) has a unique 
solution e which belongs to the class C^'^([0,T] x M^ -^ R^). 

A5 There exists a constant K > which does not depend on {t,x,y,z), 

such that ET=l\i-Jit,^,y,{^^}Zl)\{L\Piiy)\Mdy)y' < K. 

Lemma 3.5. Assume A3, A4, and A5 hold. Then the function f{t,x,y, ■ ) o 
9^^\t,x), where 6^-^\t,x) is given by (3.3), is Lipschitz in x for all {t,x,y), 
and the Lipschitz constant does not depend on {t,x,y). 

Proof. Note that by Assertions 1 and 2 of Lemma 3.1, the function c^ = 
{Ci}°li is Lipschitz in two spatial variables as an ^2-valued function. By A4, 
6 and dkO are Lipschitz. Therefore, the last summand in (3.3) is Lipschitz 
in X, and moreover, its Lipschitz constant does not depend on t by bound- 
edness of the both multipliers. Next, we already proved that the map (3.12) 
is Lipschitz. This and the boundedness of the gradient of 6 imply that the 
map 

$ : M^ ^ M^ X £2, a; f^ / ^(t, x + 6{t, x, 9{t, x),y)) - e{t, x))pi{y) v{dy) 

is also Lipschitz. Argument (3.13) implies that the Lipschitz constant of $ 

has the form K I J^\pi{y)\'^i/{dy)] where .R' is a constant that does not 
depend on i. Now A5 implies the statement of the lemma. D 

Proposition 3.6. Assume A3, A4, and A5. Then, the SDE 

At = /(s,A„^(s,A,),^W(s,A,_))rfs + ^~,a,(s,X,_,^(s,A,_))rfi7«, 



(3.21) 

where 6 is the solution to (3.2) and 6^^^ is defined by (3.3), has a pathwise 
unique cddldg adapted solution. 

Proof. We will show that 

^(A)i = x+ / /(s,A„^(s,A,),0«(s,A,_))ds+ f a{s,Xs.,e{s,X.,.))dHs 
Jo Jo 
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Is 



is a contraction map in the Banach space S with the norm ||$| 
i^supjgro^T] |$tP- Take two points Xg and X[. from 5*. For simpUcity of no- 
tation, we define 0"^ = a{s, Xs_, 9{s, Xs_)) and a'^ = a{s, X'^_, 9{s, X'^_)) . To 
estimate the difference of the stochastic integrals with the integrands a^ and 
a'g with respect to the || • ||5-norm, we apply the Burkholder-Davis-Gundy 
inequality to the martingale JQicTs — crs)dHs. We obtain that there exists a 
constant C > such that 



E sup 

re[0,i] 



(as - a',)dHs 



^CE 



(a, - cr'JdHs 



J t 



Jo * ,,^^ Jo 



CE I \\(ys-ors\\lpxe2d'S 



where [ ■ ]t and {■)t are the quadratic variation and the predictable quadratic 
variation, respectively. Moreover, we applied the identity {Hi, Hj)s = SijS and 
the decomposition [M]t = {M)t + Ut for the quadratic variation of a square in- 
tegrable martingale (i.e. a martingale Mt with sup^ jM^p < oo) in the sum of 
the predictable quadratic variation and a uniformly integrable martingale Ut 
starting at zero. Next, we note that the functions x \-^ f{s, x, 9(t, x), 9^{t, x)) 
and X I— 7- a(t,x,9{t,x)) are Lipschitz whose Lipschitz constants do not de- 
pend on t. This and the above stochastic integral estimate imply that there 
exist a constant K > such that 



E sup |^(X), - ^(X').l^ ^ KE f \X, - X'fds ^ KE [ sup |X, - X'fds. 
sefo,il Jo Jo refO.sl 



^ — j— -E' sup \Xs 

^- se[o,t] 



X 



/|2 



Iterating this n — 1 times we obtain: 

E sup |^"(X), - ^"(x); 

sG[0,i] 

Choosing n sufiicienty large so that 
traction, and thus, \l/ is a contraction as well. By the Banach fixed point 
theorem, the map \l/ has a unique fixed point in the space S. Clearly, this 
fixed point is a unique solution to (3.21). Setting X^^^ = x, and then, suces- 
sively, X^"'^ = "^(X^^'^^), we can choose cadlag modifications for each X^'^K 






< 1, we obtain that \E'" is a con- 
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Since the X'^^^'s converge to the solution X in the norm of S, X will be also 
cadlag. This cadlag solution is unique in the space S, and therefore, pathwise 
unique. D 

Introduce the space S of J-'rpredictable MP x £2-valued stochastic processes 
with the norm ||$||| = E J^ ||$s||^Qx^ ^^- Now we formulate our main result. 

Theorem 3.7. Suppose A3, A4, and A5 hold. Let Xt be the cadlag adapted 
solution to (3.21). Then, the triple {Xt,Yt, Zt), where Yt = 6{t,Xt), Zt = 
9^^\t,Xt-) with 9^^^ given by (3.3), is a solution to FBSDEs (3.1). Moreover, 
the pair of cadlag solution processes [Xf, Yt) is pathwise unique. The solution 
process Zt is unique in the space S. 

Proof. It suffices to prove that the triple (X^, Fj, Zt) defined in the statement 
of the theorem verifies the BSDE in (3.1). Application of Ito's formula to 
6{s,Xs) gives: 

e{T,Xt)-e{t,Xt)= f dse{s,x,_)ds+ [ dke{s,x,_)dx^ 

Jt Jt 

+^ I dle{s,x,_)d[{x^)\{x^)%+ J2 [0{s,x,)-e{s,x,_)-Ax^d,e{s,x,_)], 

(3.22) 

where X^ is the continuous part of Xg. Using the representation for Hg from 
Lemma 2.1 we obtain that 

d[{X'f, {X'^Yl = 2/3^'(s, Xg, e{s, Xg))ds, 

where (3^^ is given by (3.5). The forward SDE in (3.1), the relation AHg = 
Pi{ALs), obtained in [10], and representation (3.4) for the function 6 imply: 

oo 

AX, = ^ai(s,X,_,n_)A/f» = 6{s,Xs_,Ys_,ALs). (3.23) 

Next, one can rewrite the last term in (3.22) as 
J2 [e{s,Xs-+S{s,X,.,Y,^,ALs))-e{s,X,^) 

t<s^T 

-5\s,Xs-,Ys.,AL,)dkeis,Xs^ 
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where 5^ is the fcth component of 5. Define the random function 

h{s,y)=9{s,Xs.+5{s,X,_,d{s,X,,),y))-9{s,X,_) 

-6\s,X,^,e{s,X,^),y)dke{s,Xs-). (3.24) 

Note that for each fixed s G [0,T] and cj G fi, the function h satisfies condi- 
tion (2.1). Indeed, the mean value theorem, e.g. in the integral form, can be 
applied to the difference of the first two terms in (3.24). By boundedness of 
the partial derivatives dk9, it suffices to verify that 

E I \5{s,Xs.,e{s,X,_),y)\^ dsu{dy) < oo. 
The latter holds by Assumption AO and formula (3.7). Now Lemma 2 implies: 
J2 his,ALs) 



t<ss£r 

oo „y 






e{s,Xs-+S{s,X,_,e{s,Xs^),y))-9is,X, 



-S^{s,X,.,e{s,X,.),y)dk9is,X,.) pi{y)u{dy)dH^ 



(0 



+ / / 9{s,Xs- + 5{s,Xs-,d{s,Xs-),y))-9{s,Xs-) 

-6\s,Xs-,e{s,Xs-),y)dke{s,Xs-)]iy{dy)ds. 

Substituting this into (3.22), replacing dX^ with the right-hand side of (3.21), 
and taking into acccount that Yt = 6(t,Xt) and that 6{T,Xt) = h{XT) by 
(3.2), we obtain: 

Yt = h{XT)- j 9,0(s,X,_) + 9fc^(s,X,_)/'=(s,X,_,0(s,X,_),0«(s,X,_)) 

+ \due{s,X,_)fi^\s,X,_,9{s,X,_)) + [ [e{s,X,_ + 6{s,X,_,e{s,X,^),y)) 
^ Jr 

-e{s,X,^)-S\s,Xs^,e{s,Xs-),y)dke{s,Xs^)]u{dy) ds 

- I X^[/[^(s,^s- + '5(s,X,_,^(s,X,_),y))-e(s,X,_) 
Jt ^^-^ \_Jr 

- ^ke{s,Xs.)c\{Xs-.y)\v^{y>{dy)] dH^\ 
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Clearly, in the first three summands under the rfs-integral sign one can equiva- 
lently write Xg or Xs_. This is true since Xg has cadlag paths, and therefore, 
Xg and Xs_ can differ only at a countable number of points. Now taking 
into account PIDE (3.2), we note that the integrand in the drift term is 
—g{s,Xs-,9{s,Xs_),9^^\s,Xs~)) which is —g{s,Xs-,Ys-,Zs) by the defini- 
tions of Yg and Zg, or, it can be replaced by —g{s, Xg, Ys, Zg) since Xg and Yg 
have cadlag paths. Finally, by (3.3) and the definition of Zg, the integrand 
in the stochastic term is Zg. Consequently, 



Yt = h{XT)+ g{s,Xg,Yg,Zg)- ZgdK 



'S) 



which implies that {Xs,Yg, Zg) is a solution. 

Let us prove the uniqueness. Let {Xg.Yg.Zg) be an arbitrary solution 
to (3.1). Let Yg = 9(s,Xs), and Zg = 9^^\s,Xg_), where 9 is the solution to 
(3.2), and 9^^^ is defined by (3.3). By the above argument, {Xg, Yg, Zg) verifies 
the BSDE in (3.1). Applying Ito's product formula to \Yt — Yt\'^ and taking 
into consideration that Yt = Yt, we obtain: 

\Yt -Yt\^ = -2 j (y;_ - n_, d{Yg - Yg)) + [Y- Y]t -[Y- Y]t. 
Taking the expectations in the above relation gives: 

2 r ~ 

E\Yg — Yg\ +E \\Zg — ZgW^Q^^^ds 

= 2E j {Yg - Yg, g{s, X„ Yg, Zg) - g{s, Xg, Yg, Zg)) ds. 
By A3, there exists a constant C > such that 

E\Yt-Yt\' + E WZg-ZgWla^e^ds 

^CE [ \Yg-Yg\{\Yg-Yg\ + \\Zg-Zg\\^Q,,^)ds. 



Now using the standard estimates and applying Gronwall's inequality, we 
obtain that E\Yt — Yt]"^ + cE f^ \\Zg — ZgW^g^^ ds = for some constant 

c > 0. The latter relation holds for all t G [0,T]. This proves that Yt is a 
modification of Yt and that \\Z — Z\\s = 0. The uniqueness result follows. D 
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4 Option Pricing with a Large Investor in 
Levy-type Markets 

Usually, when modeling financial assets it is assumed that all investors are 
price takers whose individual buy and sell decisions do not influence the price 
of assets. Cvitanic and Ma [3] have already developed a model for hedging 
options in the presence of a large investor in a Brownian market. However, 
observation of real data suggests that patterns, like skewness, kurtosis, or the 
occurence of jumps are sufficiently significant (see, e.g., Eberlein and Keller 
[4]) to deserve to be accounted in a realistic model of option pricing. Further- 
more, the graphs of the evolution of stock prices at different time-scales are 
sufficiently different from the self-similarity of a Brownian motion. Thus, we 
develop a Levy-FBSDE option pricing model. We believe that such a model 
conveys a much more realistic approach to option pricing in the presence of 
the already mentioned empirical market characteristics. We assume the exis- 
tence of a Large investor, whose wealth and strategy may induce distortions 
of the price process. 

Let A^ be a Levy-type Market, i.e. a market where the stock price dy- 
namics St obeys the equation St = 5*06^*, where Xt is a Levy-type stochastic 
integral [1]. The market consists of d risky assets and a money market ac- 
count. For the price process Po{t) of the money market account, we assume 
that its evolution is given by the following equation 

dPo{t) = Po{t) r{t, W{t), Z{t)) dt, O^t^T, 
Po{0) = 1, 

where W{t) is the wealth process, and Z{t) is a portfolio-related process in a 
way that will be explained later. For the risky assets, we add the stochastic 
component represented by the volatility matrix a taking values in M'^ x £2- We 
postulate that the evolution of the d-dimensional risky assets price process 
P(t) = {Pi(t)}ti is given by the following SDK: 

00 
dP,{t) = fi{t, P(t), W{t), Z{t)) dt + Y, ^% Pit). Wit)) dHi'\ 

i=i 
PiiO)=Pi, p^^0, l^t^d, te[0,T]. (4.1) 

We derive the BSDE for the wealth process as in [3]. For the convinience of 
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the reader we repeat this derivation: 

dw{t) = Y^a.{t)dpm + ^^'^~^h^'^'^^'^'^ dm), 

where ai{t) is the portfolio process. Substituting dPi{t) with the right-hand 
sides of (4.1), we obtain: 

d oo 

dw{t) = J2 «^(^) [Mt^ nt). w{t), z{t))dt+Y, ^% Pit), w{t)) dH^^\t)\ 
i=i i=i 

d 

+ {Wit) - Y, ^^P^it)) rit, Wit), Zit)) dt 

oo 

= g{t, P{t), W{t), Z{t), a{t))dt + J2 Z,{t)dH^'\t), (4.2) 
where 

d d 

g{t, TT, w, z,a) = 'Y aifi{t, n, w;, z) + (w - ^ ai^i) r{t, w, z), 

a = {ai}f=i, vr = {'^iYi=i', 
d 

Z,{t) = Y,o^^{tyi{t,P{t),W{t)), 1 = 1,2,.... (4.3) 

i=i 
As we are assuming the absence of risk for the money market account, the 
evolution of its price depends totaly on the interest rate the investor is earn- 
ing. On the other hand, to describe the evolution of the risky assets we use 
an SDE with the stochastic term given as a sum of stochastic integrals with 
respect to H^^'^^s. This adds explanative power to the model, as it affords 
the isolation of the individual contributions of each H^^\ Now, to guaran- 
tee that the stock price has positive components we will rewrite (4.1) for 
Qi{t) = logPj(t) using Ito's formula. For simplicity of notation, we will use 
the same symbols /, g, a, and h for the coefficients of the FBSDEs which 
we obtain after rewriting SDE (4.1) with respect to Q(t) = {Qiit)}f^i and 
substituting Pi(t) = exp{(5i(t)}: 

Q{t) = q+ [ fis,Qis),Wis),Zis))ds+ [ ais,Qis),Wis))dHis), 
Jo Jo 

(4.4) 
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where q = {\ogpi}f^i. Due to relation (4.3), we exclude the dependence on 
a{t) in (4.2). BSDE (4.2) takes the form: 

W{t) = hiQiT))+ [ g{s,Qis),W{s),Z{s))ds- [ Z{s)dH{s). (4.5) 



Theorem 3.7 and relation (4.3) imply the following result. 

Theorem 4.1. Assume A3, A4 and A5. Then, FBSDEs (44-4.5) has a 
unique solution {Q{t),W{t), Z{t)) such that the pair {Q{t),W{t)) is cddldg. 
Furthermore, if for some M."^ -valued stochastic process {c(j{t)}'^^i, aj{t) ^ 0, 
relation (4.3) holds in the space S, then {aj(t)}'^^i is a replicating portfolio. 

It is evident that Levy-type markets pose theoretical questions that had 
never been raised in the Brownian motion framework. We conclude by reen- 
forcing the idea that the impossibility of replicating every potential contin- 
gent claim is rather an expected characteristic due to the complexity of the 
price formation in Levy-type markets, than a drawback of our model. 
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